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Abstract. In my review of jet phenomena in quasars, I focus on the 
following questions: How powerful are jets in radio-loud quasars? What 
is their composition? How are they launched? And why, in most quasars, 
are they so weak? I demonstrate the exceptional role that blazar studies 
can play in exploring the physics and structure of the innermost parts of 
quasar jets. 



1. Introduction 

The jet activity in quasars is common, but very diverse. As radio observations 
indicate, jet powers can differ by several orders of the magnitude within the same 
optical luminosity range. The most powerful jets produce hundred kiloparsec- 
scale double radio structures. They can be characterized as composed from a 
pair of edge-brightened radio lobes, with the hotspots matching their luminosity 
peaks. Additionally, one sided jets are often observed, connecting one of two 
hotspots with the radio core in the center of the host galaxy. The above struc- 
ture has a good interpretation in terms of a dynamical model which involves 
propagation of light, relativistic jets in the IGM (Scheuer 1974; Begelman & 
Com 1989). The hotspots are located at the ends of channels drilled by the jets 
through the IGM. They mark the regions where material of the jet is shocked 
and spreads sideways, forming the radio lobe. Relativistic speeds of jets explain 
their one-sided appearance, while lightness (jet density lower than IGM den- 
sity) is necessary to explain formation of extended radio lobe structures and the 
non-relativistic speeds of hotspots. 

Highly polarized and relatively steep radio spectra of extended radio struc- 
tures are uniquely interpreted in terms of optically thin synchrotron radiation. 
The synchrotron spectra extend from ~ 10 MHz, up to IR/optical. I discuss 
briefly, in §2, how these data can be used to estimate the jet power. 

Quasar jets can be traced in radio down to parsec-scale central regions. 
There the jets are more relativistic (bulk Lorentz factor T ~ 10 [Padovani & 
Urry 1992; Ghisellini et al. 1993; Homan et al. 2000]) than on kpc-scales (r ~ 3 
[Wardle & Aaron 1997]), and only those which are oriented close to the line of 
sight are bright enough to be observed in detail. Parsec-scale jets viewed "pole- 
on" can be decomposed into radio-cores and one-sided linear structure. The 
linear structure is very inhomogeneous, with some bright regions propagating 
with relativistic speeds and appearing to us as "superluminal" sources. Parsec- 
scale radio sources show spectra with a low-energy break due to synchrotron- 
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self-absorption. The frequency of this break is larger the closer to the center 
one measures, and superposition of spectra from all radio components give the 
characteristic flat radio spectrum, with a energy spectral index a < 0.5. Quasars 
with such radio spectra are called FSRQ (flat-spectrum-radio-quasars). I discuss 
the energetics and composition of quasar parsec-scale jets in §3. 

As high frequency VLBI observations of nearby radio galaxies show, extra- 
galactic jets are launched much deeper than the angular-resolution and synchrotron- 
self-absorption limited observations can follow in quasars (Lobanov 1998; Junor, 
Biretta, & Livio 1999). Fortunately, sub-parsec scale jets radiate a lot at higher, 
non-radio frequencies: up to optical/UV by the synchrotron mechanism, and in 
the X-ray and 7-ray bands via Comptonization of synchrotron and external dif- 
fuse radiation fields (Sikora, Begelman, & Rees 1994; Blazejowski et al. 2000). 
This radiation, Doppler boosted into our direction, often dominates over thermal 
quasar components, such as: UV/optical radiation of the accretion disk, X-ray 
radiation of a disk-corona, and IR radiation of dust located in a molecular torus 
and heated by a disk. FSRQ which have spectra dominated by the non-thermal 
radiation from a jet are called blazars. Historically, this category also includes 
BL Lac objects, those sources with thermal signatures too weak present at all 
to place them in the quasar category. Hence, in order to avoid confusion while 
talking about "quasar-hosted" blazars, I will call them "Q-blazars". 

The broad-band spectra of Q-blazars can in general be superposed from the 
radiation produced over a large distance range. However, at least during short- 
term high amplitude flares the spectra are dominated by radiation produced 
co-spatially, very likely in short-lived shocks, somewehere at 0.1-1.0 pc from the 
center. Thus, as discussed in §4, multiwavelength studies of flares in blazars 
provide exceptional tools for exploring the structure and physics of jets on sub- 
parsec scales. 

Jets are predicted to produce radiation not only by relativistic electrons 
(here and after the term 'electrons' is used for both electrons and positrons), 
but also by cold electrons. Streaming with a bulk Lorentz factor T ~ 10, the 
cold electrons Compton scatter the external optical/UV photons and boost them 
up to the soft X-ray range (Begelman & Sikora 1987; Sikora et al. 1997). It 
should be emphasized that cold electrons are an unavoidable constituent of jets 
near their base, where even mildly relativistic electrons cool faster than they 
propagate. And they can be present up to distances where non-thermal flares 
are produced, dragged by as-yet-unshocked portions of the flow. As yet, no 
soft X-ray excesses have been confirmed. The upper limits imposed on the 
number of cold electrons by observed soft X-ray fluxes exclude pure e + e~-jets 
and provide strong constraints on the minimum distance of jet acceleration and 
collimation (Sikora h Madejski 2000). These constraints, together with possible 
jet production scenarios, are discussed in §5. 

Of course, any model of jet production should be able to explain the huge 
range of jet powers. Recent discoveries that many luminous radio-quiet quasars 
reside — like the radio-loud quasars — in giant ellipticals (Taylor et al. 1996; 
Kukula et al. 2000), and that the galaxy environments of the same luminosity 
radio-quiet and radio- loud quasars is similar (McLure & Dunlop 2000) , challenge 
the previous claims that radio-loudness can be related to the morphology of 
the host galaxy or its clustering richness. Furthermore, optical/UV spectral 
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similarities (Francis, Hooper, h Impey 1993; Zheng et al. 1997) and recent 
discoveries that BAL (broad absorption line) systems exist also in radio-loud 
quasars (Brotherton et al. 1998; Becker et al. 2000) suggest that radio-loudness 
is not very dependent on the parsec-scale environment, as well. All the above 
strongly supports the so-called spin paradigm, according to which powerful jets, 
giving rise to radio-loud quasars, can be produced only with the help of rapidly 
rotating black holes. How the spin paradigm relates to central engine models 
and the evolution of quasars is discussed toward the end of §5. 

2. Radio Lobes 

What can we learn about jets from radio lobes? First of all, they provide very 
useful information about the energetics of jets, which, contrary to that derived 
from radio-core scales and smaller, is not biased by such uncertainties as jet 
bulk Lorentz factor, dissipation efficiency and variability. The procedure for de- 
riving the jet power is simple, but not free of assumptions and approximations. 
The first step is to recover from the observed electromagnetic spectrum the en- 
ergy distribution of electrons. This can be done using the following approximate 
formula: 

Lu,syndv ~ (N^d^meC 2 ^ (1) 

where 

4 ca T u B l 2 rn , 

7 - o — 2 

3 m e c z 

is the rate of electron synchrotron energy losses; ub = B 2 /8it is the magnetic 
energy density; and u oc ^ 2 B. For power-law synchrotron spectrum L u cx v~ a , 
formula (1) gives iV 7 oc j~ s , where s = 2a + 1. The above procedure is not 
sufficient to determine the normalization of the electron energy distribution, 
however. Additionally, one needs to know the intensity of the magnetic field, 
which can be estimated by assuming energy equipartition between electrons 
and magnetic fields. With this assumption, the electron (and magnetic) energy 
content of the quasar radio lobes is found to be in the range 10 59 — 10 61 ergs 
which, when divided by spectrally or dynamically determined ages of the radio 
lobes, tiobe ~ 3 x 10 7 years, gives jet powers 10 45 — 10 47 ergs s _1 (Rawlings and 
Saunders 1991). 

There are several reasons why the above estimates should be considered 
as lower limits. First, the equipartition condition corresponds almost exactly 
with the minimum total energy of electrons and magnetic fields. Second, from 
the observed radio spectra one can deduce that most of the energy carried by 
electrons is contained in the low energy part of their distribution. Assuming 
that the electron distribution has a break at an energy corresponding with the 
lowest observable frequency ~ 10 MHz (limited by reflection of radio waves 
in ionosphere), one finds that ^ m in ~ 500/^B/lOfiG. Since there is no proof 
that radiation has an intrinsic cutoff at 10 MHz, the adopted value of j m i n can 
be greatly overestimated and the total electron energy content underestimated. 
Third, Rawlings & Saunders assumed no energy contribution from protons. 

How much might the jet powers be underestimated due to the above as- 
sumptions? Detection of X-rays from hotspots of several nearby radio galaxies 
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and their interpretation in terms of the SSC (synchrotron-self-Compton) pro- 
cess allowed one to derive the magnetic field and electron energy densities with- 
out assuming equipartition (see Wilson, Young & Shopbell 2000 and references 
therein). Departure from the equipartition condition have been found to be 
very small. However, recent observations of X-rays produced around quasar 
radio lobes show that the pressure of the external gas is several times larger 
than the pressure in the radio lobes obtained assuming equipartition between 
electrons and magnetic fields and no proton contribution (Hardcastle & Worral 
2000). This inconsistency cannot be resolved, even assuming that the electron 
distribution extends down to j m i n ~ 1, if the equipartition condition is kept. 
There must be significant departure from equipartition between electrons and 
magnetic fields and/or the lobe pressure is dominated by protons. There are 
some observations suggesting that, indeed, the equipartition conditions can be 
violated in radio lobes, with particle pressure dominant over magnetic pressure 
(see Blundell and Rawlings 2000 and references therein). Yet another possibil- 
ity is that extra pressure in radio lobes is provided by cosmic rays accelerated 
via the Fermi process operating in the boundary layer between a jet and the 
surrounding medium (Ostrowski 2000). 

The rate at which energy is delivered to radio lobes can be estimated more 
directly, just from the bolometric luminosities of hotspots. The rate is 

Lj = ^1. (3 ) 

where r] e is the fraction of kinetic energy of a jet converted in the shock to 
relativistic electrons and rj ra d is the fraction of electron energy lost by radiation. 
If, in the radiative regime, the electromagnetic spectrum has a slope a ~ 1, 
then ri rad ~ \n(v max / v c ) / ln(v max / v min ) , where v c is the "cooling" break, and 
v-min < 10 MHz. Applying this for hotspot D in Cyg A, where v c ~ 10 10 
Hz and v max ~ 10 12 Hz, one can find that for v m in = 10 MHz, rj ra d ~ 0.4. 
Now, assuming that the energy dissipated in the hotspots is equally shared by 
electrons and magnetic fields, i.e., rj e = 1/2, and adopting from Meisenheimer 
et al. (1997) L HS ~ 4 x 10 44 ergs s _1 , we obtain L j ~ 2 x 10 45 ergs s _1 . This 
estimate of Lj is consistent with that deduced from the radio-lobe energetics 
(Carilli Sz Barthel 1996). In both cases the energy is underestimated only by a 
factor 3/2, if there are protons and they equally share energy with electrons and 
magnetic fields. Unfortunately, hotspot spectra up to the highest synchrotron 
frequencies are currently available only for nearby radio galaxies and, therefore, 
the bolometric-luminosity method cannot be applied to distant quasars. 

What about the pair content of radio-lobe plasmas? Noting that jets are 
approaching the hotspots with relativistic speeds, the average energy of shocked 
protons is expected to be of the order of the jet Lorentz factor. Thus, if the dis- 
sipated kinetic energy of a jet is initially shared equally by electrons, protons, and 
magnetic fields, the average electron energy should be 7 ~ (l/3)T(m p /m e )(n p /n e ) 
~ 600r (n p /n e ). Since radio observations cannot follow electrons with energies 
lower than 7 ~ 500/ ^(B /lOfiG), the pair content cannot be verified by radio 
lobe observations. 
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3. Parsec-scale Jets 

In order to recover the physical parameters of radiating plasma in compact radio 
sources, one needs: to use the emissivity formula (like that in equation [1], but 
written in the source comoving frame); to transform the comoving luminosity 
and frequency to the observed ones; and to take advantage from two of the 
following: 

a) the value of the bulk Lorentz factor T (if available from VLBI observa- 
tions); 

b) X-ray flux, provided the X-rays are produced by the SSC process; 

c) synchrotron-self-absorption break; 

d) equipartition condition. 

Such analyses have been performed for a large sample of compact radio sources 
by Ghisellini et al. (1992), who used (b) and (c), and for the series of compact 
radio components in 3C 345 and 3C 279 by Hirotani et al. (1999; 2000), who used 
(c) and (d). Among other aspects, they calculated electron energy distributions, 
n 7 = Cj~ s , assuming ~f m in = 1> where s = 2a + 1. This allows one to obtain 
electron energy densities, u' e = v!^)m e c 2 , and then energy fluxes of relativistic 
electrons 

L e ~ u' e r 2 Tra 2 c (4) 

where a is the cross-sectional radius of the source. For the studied sources, 
values of L e are in the range 10 45 — 10 47 ergs s -1 , provided j m i n = 1, and 2-3 
times smaller, if r y m in 3> 1. 

Hence, the energy fluxes of relativistic electrons alone come close to satisfy- 
ing the energy requirements of radio lobes. However, one should note that the 
total energy flux also includes other forms of energy and, in general, we have: 

Lj = ^ r (5) 

WdissVeil ~ Vrad) 

where rjdi ss is the fraction of the total jet energy which is dissipated and used 
to accelerate electrons, to heat protons and to amplify magnetic fields; r} e is the 
fraction of the dissipated energy which is used to accelerate electrons; and r) ra d 
is the fraction of electron energy lost by radiation during the source lifetime. 
Thus, with the derived values of L e , the total energy fluxes, Lj, become danger- 
ously high, particularly if rjdiss < 0.1 as intrinsic shock theories predict. There- 
fore, provided that the derived densities of relativistic electrons are not affected 
by systematic errors (noting their very strong dependence on the absorption- 
turnover frequency and the source geometry), one needs to postulate external 
shock models, with dissipation efficiencies r]di ss > 0.5 (Dermer & Chiang 1998). 

What is the pair content of the compact radio sources? Assuming that the 
dominant energy carriers are protons, we have 

Lj ~ L P: q + Ldiss = L Pi q + 5L P + Lb + L e , (6) 

where L Pi o = n' p m p c 3 TTa 2 T 2 and 5L P ~ n' p m p {i^ p ) — l)c 3 r 2 7ra 2 . Since 

L e n' e {~i)m e 7] e 7] dis s(l - ri rad ) 



Lpfi n' p m p 1 - rjdi 



(7) 
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where r\ e = 1/(1 + Lb/ L e + 5L p /L e ), we obtain 

4 - ^ - — (3^-^Ml - rw) (8) 

Tip ZTlp ^fmin J- Vdiss 

where I used (7) ~ 6j m i n , which corresponds with a ~ 0.6. 

Unfortunately, due to the synchrotron-self-absorption it is impossible to 
follow electrons with energies 7 < 50^ (f a 5 s /lGHz) /(£?/0.1G) and determine 
^min directly from the observed spectra. One can eventually try to determine 
the upper limits for the minimum electron energies from measurements of circu- 
lar polarization provided such polarization results from the Faraday conversion 
mechanism. Measurements of circular polarization have been completed for sev- 
eral extragalactic sources (Homan & Wardle 1999) and at least in 3C 279 there 
are indications that circular polarization is produced by this mechanism and 
that -y min < 20 (Wardle et al. 1998). For such 7mm 's the number of pairs per 
proton can range from a few up to tens. There are several indirect arguments in 
favor of rather low pair contents: (i) if pairs are produced in the central engine or 
its vicinity, then their flux is very limited by the annihilation process (Ghisellini 
et al. 1992); (ii) if pairs are created by nonthermal pair cascades operating in 
the jet shocks, they would produce much softer X-ray spectra than observed in 
Q-blazars (Ghisellini & Madau 1996); (iii) jets with a large number of cold elec- 
trons would produce soft X-ray bumps by Comptonization of UV disk radiation 
and BELs and such bumps have not been confirmed (Begelman & Sikora 1987; 
Sikora et al. 1997). 



4. Q-blazars 

Strong and fast variability in blazars is commonly interpreted in terms of the 
shock-in-jet model. Producing a flare of the observed time scale tfi, the shock 
passes a distance range 

Ar fi ~ ct fl T 2 ~ 100(t // /lday)(r/10) 2 It - days. (9) 

Sharp profiles of flares and comparable time scales of their rise and decay suggest 
that shocks are launched at distances 777 ~ A777 (Sikora et al. 2000). The 
distance of flare production can be also estimated from the spectral location 
of the 7-ray luminosity peak, provided the peak is related to the break in the 
electron energy distribution caused by the cooling effect, i.e. that above the 
peak electrons radiate on time scales shorter than the lifetime of the shock. If 7- 
ray production is dominated by the ERC (external-radiation-Compton) process, 
and the luminosity of the ERC component is larger than the luminosity of the 
synchrotron component, then the spectral distance is 

r sp ~ ct r ~ J —U- (10) 

v o T V v c u diff 

where the following relations were used: t' = |7 C /7 C |; [7] ~ T 2 ^ 2 UTUdif / / 'm e c; 
and 7 C ~ J u c /u,nff /T . It can be checked that both the variability distance, rfi, 
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and the spectral distance, r sp , are of the same order if Udiff ~ 0.005 ergs cm 3 , 
i.e., if LsELiffi) ~ 10 44 ergs s _1 , and/or covering of the central source by dust 
at T = 1000 K is ~ 0.1. Both quantities are consistent with ourknowledge about 
BELR and near-IR radiation in quasars (Blazejowski et al. 2000). 

Equations (9) and (10), combined with emissivity formulae for the produc- 
tion of X-rays via the ERC process (Blazejowski et al. 2000) and with given 
value of Ldiff(rfi), can be used to calculate the number of relativistic electrons 
involved in flare production and their energy flux, L e . For Lbel ~ 10 45 ergs s _1 
and T calculated from the model, it can be found that L e ~ 10 L$x,46 ergss -1 , 
where L$x is the soft X-ray luminosity (Sikora et al., in preparation). Since 
soft X-rays are very likely dominated by the SSC process (Inoue & Takahara 
1996; Kubo et al. 1998; Blazejowski et al. 2000), the above estimate should be 
considered only as the upper limit. 

The electron energy flux, L e , can be estimated also using the bolometric 
luminosity procedure. We have 

L e ~-^ (11) 

4vr Tjrad 

where Q, ~ ir/T 2 and Lqb is the total apparent luminosity of a blazar. During 
high states Lqb is dominated by luminosity in 7-ray bands and is of the order 
10 48 ~ 49 ergs s -1 (von Montigny et al. 1995). From typical high energy spectra of 
Q-blazars (a ~ 1 for hu > 30MeV and location of v c in the 1 — 30 MeV range), 
one can conclude that rj ra d ~ 1/2 if 7 m j n ~ 1, and rj rac i ~ 1 if 7 m j„ ~ 100. 
With these numbers, and T taken from the model calculations, we obtain L e ~ 
1— 2xl0 45 (r/15) _2 LQs ! 48 ergs s -1 , which is of the same order as calculated using 
the X-ray flux. Noting that in Q-blazars L e is about 1-2 orders lower than Lj 
estimated from radio lobes, and that L e = rjdi ss VeLj(l — rj ra d), one can conclude 
that activity of the sub-parsec jets is governed by low dissipation efficiencies, 
very likely by intrinsic shocks. This conclusion can actually be reinforced by the 
fact that the high (flaring and radio active) states are not permanent, and after 
averaging over longer periods of time the values of L e are likely to be lower by 
at least a factor two. 

If the main carriers of energy are cold protons, then from equation (8) the 
pair content is n' + /n' p ~ 5/^ min for j min < 5, and is negligible for ^ rnin > 5, 
assuming a = 0.6, i]di ss = 0.05 and r] = 1/3. Here, like in the radio lobes and 
compact radio sources, we have a sort of conspiracy regarding the value of ^ m in- 
Since X-rays from blazars are presumably dominated by the SSC process, it 
is very difficult to follow the low energy portions of the ERC spectra to check 
whether there are any signatures of the low-energy cutoff in energy distribution 
of electrons. But at least we can say that a large pair content n+/n p 3> few 
can be excluded. One could speculate that there is a large number of cold 
pairs. However, this is excluded because such pairs, Comptonizing external 
UV photons, would produce a huge soft X-ray bump, which is not confirmed 
observationally. 
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5. Central Engine and Spin Paradigm 

The absence of soft X-ray bumps in Q-blazar spectra provides strong constraints 
on the jet's structure near its base. There, Comptonization of the disk radiation 
by cold electrons, dragged by the relativistic and well collimated jet, is predicted 
to produce a prominent X-ray bump even if number of pairs is zero. This suggests 
that powerful jets are wider and/or slower at their bases (Sikora & Madejski 
2000). They can be launched by the innermosts parts of the accretion disk, with 
the matter pulled from the disk surface and accelerated along the open magnetic 
field lines by centrifugal forces (Blandford & Payne 1982). Such proto-jets are 
probably collimated further away, by a disk corona or winds predicted by some 
models to be formed at distances > 100i? 9 (Rozanska & Czerny 2000; Murray 
& Chiang 1997; Proga, Stone & Rahman 2000). Yet in the acceleration zone, 
the initially proton-electron outflows can be loaded by e~-e + pairs. This is due 
to boosting of the coronal hard X-rays by cold electrons in the outflow up to 
MeV energies and the subsequent absorption of MeV photons in the 77 pair 
production process. Depending on the geometry and kinematics of the proto- 
jet, the number of pairs per proton can reach a value ranging from a few up to 
a few tens. 

The fact that only 10 % of quasars are radio-loud implies that powerful jets 
are rare. It is very likely that the leading parameter which decides about the jet 
power is the spin of the black hole, defined as A = J/J max where J max = GM 2 /c 
(Wilson & Colbert 1995; Moderski, Sikora & Lasota 1998). Perhaps the fast 
rotation of the black hole is necessary to heat the accretion disk. Due to this 
heating the innermost parts of the disk can be inflated and this can help to 
generate large scale magnetic fields, which are required both to accelerate MHD 
outflows and to link a disk with the rotating black hole and/or the gas plunging 
into the ergosphere (Meier 2000; Krolik 2000). Furthermore, extra heating of 
surface layers of the disk by a rotating black hole can help to put gas on magnetic 
field lines, if the latter are not bent enough to allow centrifugal forces to pull 
the gas directly from the photosphere. 

Now, provided that the above scenario(s) gives a strong enough dependence 
of the jet power on the spin of the black hole to explain the huge range of radio- 
loudness, the next question which should be answered is how Nature managed 
to have only a small fraction of rapidly rotating black holes. Let us recall that 
accretion disks, acting enough long to double the black hole mass, spin up the 
holes to A ~ 1 (Bardeen 1970; Thorne 1974). On the other hand, once the black 
hole is spun up, it cannot be easily slowed down, certainly not during the low 
accretion phases during which the rate of extraction of black hole energy is very 
small because of its proportionality to B 2 oc M. Hence, if the spin paradigm 
is right and the growth of supermassive black holes is governed by the large 
accretion events, then the number of radio loud quasars should be larger than of 
radio quiet quasars, oppositely to what is observed. A solution of the problem 
can be that growth of the black hole in most objects is dominated by low-mass 
accretion events with random angular momentum orientations. This picture 
is supported by observations of Seyfert galaxies, which show that AGNs are 
randomly oriented relative to the host galaxy planes (Wilson & Tsvetanov 1994; 
Schmitt et al. 1996). Low- mass accretion events are also supported by analyses 
of the dynamics of Seyfert jet interactions with the host galactic gas which show 
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that the lifetime of these objects is of the order of 10 5 years only (Capetti et 
al. 1999). Thus, it is tempting to speculate that only major mergers which 
involve at least a one gas-rich galaxy can lead to the accretion disk operating 
long enough to double the black hole mass and spin up the black hole to A > 0.5. 
During early phases of the process, when A is still small, the quasar is predicted 
to show up as radio-quiet, and then slowly transform to a radio- loud one. After 
the fuel is off, the accretion drops, but spin of the black hole remains high and 
such an object can eventually be represented by FR I radio galaxies. 

Acknowledgments. I am grateful to Mitch Begelman for his valuable com- 
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by the Polish KBN grant 2P03D 00415. 
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